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Abstract. Weakly interacting massive particles (WIMPs) are one of the leading candidates for Dark 
Matter. So far we can use direct Dark Matter detection to estimate the mass of halo WIMPs only by 
fitting predicted recoil spectra to future experimental data. Here we develop a model-independent 
method for determining the WIMP mass by using experimental data directly. This method is 
independent of the as yet unknown WIMP density near the Earth as well as of the WIMP-nuclear 
cross section and can be used to extract information about WIMP mass with 0(50) events. 
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1 Introduction 

Today astrophysicists have strong evidence [l]-[5] to 
believe that a large fraction (more than 80%) of the 
matter in the Universe is dark (i.e., interacts at most 
very weakly with electromagnetic radiation) . The dom- 
inant component of this cosmological Dark Matter (DM) 
must be due to some yet to be discovered, non-baryonic 
particles. Weakly interacting massive particles (WIMPs) 
X with masses between 10 GcV and a few TeV are one 
of the leading candidates for DM [6] . 

Currently, the most promising method to detect 
different WIMP candidates is the direct detection of 
the recoil energy deposited in a low-background labo- 
ratory detector by elastic scattering of ambient WIMPs 
off the target nuclei [7], [8]. The differential rate for 
elastic WIMP-nucleus scattering is given by [6] : 
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Here R is the direct detection event rate, i.e., the num- 
ber of events per unit time and unit mass of detector 
material, Q is the energy deposited in the detector, 
F{Q) is the elastic nuclear form factor, and fi{v) is 
the one-dimensional velocity distribution function of 
the WIMPs impinging on the detector. The constant 
coefficient A is defined as 
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where po is the WIMP density near the Earth and 
(To is the total cross section ignoring the form factor 
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is defined by 
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where rriy^ is the WIMP mass and ton that of the target 
nucleus. Finally, Vmin is the minimal incoming velocity 
of incident WIMPs that can deposit the energy Q in 
the detector: 

Vmin = a\/Q, (4) 

where we define 



(5) 



So far most theoretical analyses of direct WIMP 
detection have predicted the detection rate for a given 
(class of) WIMP(s), based on a specific model of the 
galactic halo. This can be used to estimate the mass of 
halo WIMPs only by fitting the predicted recoil spec- 
tra to future experimental data, e.g., [9]. The goal of 
our work is to develop methods which allow to extract 
information on halo WIMPs by using the experimental 
data directly. In our earlier work [10] we used a time- 
averaged recoil spectrum, assumed that no directional 
information exists, and derived an expression for es- 
timating moments of the normalized one-dimensional 
velocity distribution function of halo WIMPs: 
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where Qtiuc is the threshold energy of the detector, 
Tthrc = ('^^/c^Q)Q=Qthrc i^ estimated value of the 



Cosmology and Astrophysics 



Contributed Talk 



scattering spectrum at Q = Qthrc- In can be either 
determined from a given expression (e.g., a fit to data) 
for the recoil spectrum: 



Qmax = 200 keV, n = 1 



(7) 



/Q..,e FHQ) \dQ 
or estimated directly from the measured recoil energy 
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where the sum runs over all events in the data set. Note 
that all these expressions are independent of the as yet 
unknown WIMP density near the Earth as well as of 
the WIMP-nucleus cross section. (More details about 
the reconstruction of the velocity distribution function 
of halo WIMPs and the estimate of its moments as well 
as all formulae needed can be found in Ref. [10].) 



2 Determining the WIMP mass 

In this section we present a method for determining 
the WIMP mass based on the estimate of the mo- 
ments (u") from two (or more) experimental data sets 
with different target materials. The basic idea is that, 
from independent direct WIMP detection experiments 
with different target nuclei, the measured recoil spec- 
tra should lead to the same (moments of the) veloc- 
ity distribution function of incident WIMPs, when the 
threshold energies of these experiments are low enough. 
Note that this can be done independent of the detailed 
particle physics model, which determines the value of 
(To for the target nuclei. However, one will need to 
know the form factor of each target nucleus, which 
strongly depends on whether spin-dependent or spin- 
independent interaction dominates. 

Moreover, as shown in Ref. [10], because data bin- 
ning is not required for estimate of the moments (w"), 
some non-trivial information can already be extracted 
from 0(20) events, it thus seems reasonable to expect 
that one can also obtain meaningful information about 
the WIMP mass with pretty few events. 



2.1 Neglecting threshold energies 

For the case that the threshold energy Qthre can be ne- 
glected, the n— th moment of the velocity distribution 
function, (w"), in Eq.(6) can be reduced to 
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where X and Y are two target nuclei, In,x, Iq.x and 
so on can be estimated by Eq. (7) or (8) . Note that the 
form factor F'^{Q) in Eq.(7) or (8) for estimating In.x 
and are different. According to the definition of 
a in Eq.(5) with the expression of the reduced mass 
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Fig. 1. The curves show the ratios of the reproduced 
WIMP masses estimated by Eq.(lO) with different combi- 
nations of target nuclei to the input (true) one as functions 
of the input WIMP mass. The theoretical predicted recoil 
spectrum for the shifted Maxwellian velocity distribution 
function [6], [10] with the Woods-Saxon form factor [11], 
[6] has been used (rriN = 70.6 GeV/c^ for ^^Ge, vq = 220 
km/s, Ve = 231 km/s). The solid (blue) line, the dashed 
(black) line, and the dash-dotted (red) line are for ''"Ar + 
2®Si, ^^Ge + '"'Ar, and ^«Ge + ^^Si combination, respec- 
tively. 



mr,N in Eq-(3), and using some simple algebra, one can 
solve the WIMP mass as 
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Fig. 1 shows the ratios of the reproduced WIMP 
masses estimated by Eq.(lO) with different combina- 
tions of target nuclei to the input (true) one as func- 
tions of the input WIMP mass. ^^Si, ''"Ar, and ^^Ge 
have been chosen as three target nuclei and thus three 
combinations for 7?,„ defined in Eq.(ll) with n = 1 are 
shown. TZn has been estimated by the integral form 
of /„ in Eq.(7) with a maximal measuring energy of 
200 keV. In Fig. 1 we can see obviously a deviation 
of the reproduced WIMP mass from the input (true) 
one as input to^^60 GeV/c^. The heavier the nuclear 
masses of two target nuclei, the larger the deviation 
from the true WIMP mass. This is caused by intro- 
ducing the maximal measuring energy for estimating 
/„. For n — 1 and input — 200 GeV/c^, the devi- 
ation with Qmax — 200 keV is around 20%. However, 
for input ^120 GeV/c^, this deviation will be less 
than 5%, and if Q^ax = 250 keV or 300 keV, this 
deviation will be reduced to 10% or even only 5%. 
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Fig. 2. The curves show the statistical errors estimated 
by Eq.(12) with diflterent combinations of target nuclei as 
functions of the input WIMP mass. Each experiment has 25 
events, i.e., totally 50 events. Parameters and indications 
of the lines as in Fig. 1. 



Furthermore, the statistical error on the reproduced 
WIMP mass can be obtained from Eq.(lO) directly as 



Omax ^ 200 keV, Qmin ^ 1 keV, n ^ 1 , 250 + 250 events, Ge-76 + Si-28 
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The formulae for estimating tT^(/„.x) = cav{In,x,In,x) 
cov{Io^x,In,x) and so on can be found in Ref. [10]. 

Fig. 2 shows the 1— cr statistical errors estimated by 
Eq.(12) with three different combinations of target nu- 
clei as functions of the input (true) WIMP mass. Each 
experiment has 25 events, i.e., totally 50 events. Note 
that, in order to estimate cov(/„,/m) above, a thresh- 
old energy Qmin — 1 keV for both nuclei has been 
given. In Fig. 2 we can observe that the larger the mass 
difference between two detector nuclei, the smaller the 
statistical error will be. Hence, the combination with 
the largest mass difference, ^^Ge -|- ^*Si, will have the 
smallest statistical error. On the other hand, despite 
of the factor l/\n\ in Eq.(12), it has been found that 
the statistical errors increase generally with increasing 
n. Hence, n 

<T{m^) in Eqs.(lO) and (12), respectively. 

Figs. 3 show the reproduced WIMP mass with the 
1— (7 statistical error by using ^^Ge and ^^Si as two 
target nuclei as a function of the input (true) WIMP 
mass. From the upper frame, it can be found that, de- 
spite of the very few (25 -I- 25, totally 50) events and 
correspondingly very large statistical error, for < 
100 GeV/c^, one can already extract some meaning- 
ful information on the WIMP mass. For example, for 
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Fig. 3. The reproduced WIMP mass with the statistical 
error by using ^^Ge and ^*Si as two target nuclei as a func- 
tion of the input WIMP mass. The solid (red) line indi- 
cates the reproduced WIMP mass estimated by Eq.(lO), 
the dashed (red) lines indicate the 1—a statistical error es- 
timated by Eq.(12). The straight dash-dotted (green) line 
indicates the input (true) WIMP mass. Each experiment 
has 25 (250) events, i.e., totally 50 (500) events, in the 
upper (lower) frame. Parameters as in Fig. 1. 



= 25 GeV/c^ and m^^ = 50 GeV/c^, we will re- 
produce ruy ~ (25 ± 13) GeV/c^ and ~ (50 ± 
31) GeV/c\ For the case with 500 (250 + 250) total 
events, the statistical error will be reduced to less than 
5 and 10 GeV/c^, respectively! Certainly, as shown 
in the lower frame of Figs. 3, for the case with 500 
total events, the deviation of the reproduced WIMP 
mass from the input one becomes important. Never- 
theless, in practice, an experiment with more than 200 
events should have a larger maximal measuring energy, 
and, as discussed above, the deviation can (should) be 
strongly reduced. 

For the simplified simulations with the integral form 
of /„ presented above, the event numbers from both 
experiments have been considered to be equal. Practi- 
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cally experiments with the higher mass mielei, e.g.. Ge 
or Xe, are expected to measured (much) more signal 
events. However, some detailed simulations show that, 
the reproduced WIMP mass with the statistical un- 
certainty estimated by two experiments with different 
maximal measuring energies and different exposures, 
or, equivalently, different event numbers, will be mod- 
ified only slightly from that shown in Figs. 3. More- 
over, a detailed analysis of contributions from different 
terms of cr(m^) show that one can not reduce the sta- 
tistical error of estimated by Eq.(12) by improv- 
ing only one experiment with even very large event 
number, since the contribution from the other (poor) 
experiment will dominate the error. 



2.2 Taking into account threshold energies 

For the case that Qthre in Eq.(6) can not be neglected, 
TZn defined in Eq.(ll) should be modified to the fol- 
lowing general form: 



affect the reproduced WIMP mass very significantly. 
Moreover, the numerical analysis shows also that, for 
WIMP masses < 100 GeV/c^ some meaningful infor- 
mation on the WIMP mass can already be extracted 
from 0(50) total (each experiment 0(25)) events. 

The analyses of this work are based on several sim- 
plified assumptions. First, all experimental systematic 
uncertainties, as well as the uncertainty on the mea- 
surement of the recoil energy have been ignored. Com- 
paring with large statistical uncertainty this should 
be a quite good approximation. Second, the analysis 
treats each recorded event as signal, i.e., background 
has been ignored altogether. This may in fact not be 
unrealistic for modern detectors. Third, each detector 
consists of a single isotope. This is quite realistic for 
the semiconductor detectors. For detectors containing 
more than one nucleus, by simultaneously measuring 
two signals, one might be able to tell on an event-by- 
event basis which kind of nucleus has been struck. 
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where n ^ 0. In this general form of TZn there are 
totally six variables: rthrcX, In.x, lo.x and the other 
three for nucleus Y. This should generally produce a 
larger statistical error than that estimated by Eq.(12) 
due to the contribution from rthro- However, one can 
practically reduce the number of variables by choosing 
n = -1: 
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3 Summary and Conclusions 

In this paper we have presented a method which allows 
to extract information on the WIMP mass from the 
recoil energy measured in elastic WIMP-nucleus scat- 
tering experiments directly. In the long term this in- 
formation can be used to constrain e.g., SUSY models 
in the elementary particle physics and compare with 
information from future collider experiments. 

Our method for determining the WIMP mass by 
combining two (or more) experiments with different 
detector materials is independent of the as yet un- 
known WIMP density near the Earth as well as of 
the WIMP-nucleus cross section. The only informa- 
tion which one needs is the measured recoil energy. 

Due to the maximal measuring energy of the detec- 
tor, there will be a deviation of the reproduced WIMP 
mass from the true one. Nevertheless, for experiments 
with very few events and thus a quite large statisti- 
cal error in the near future, this deviation should not 



References 

1. F. Zwicky, Helv. Phys. Acta 6, (1933) 110; S. Smith, 
Astrophys. J. 83, (1936) 23. 

2. V. C. Rubin and W. K. Ford, Astrophys. J. 159, (1970) 
379; S. M. Fabcr and J. S. Gallagher, Annu. Rev. As- 
tron. Astrophys. 17, (1979) 135; V. C. Rubin, W. K. 
Ford, and N. Thonnard, Astrophys. J. 238, (1980) 471; 
K. G. Begeman, A. H. Broeils, and R. H. Sanders, Men. 
Not. R. Astron. See. 249, (1991) 523; R. P. Oiling 
and M. R. Merrifield, Mon. Not. R. Astron. Soc. 311, 
(2000) 361. 

3. M. Fich and S. Tremaine, Annu. Rev. Astron. Astro- 
phys. 29, (1991) 409. 

4. WMAP GoUab., D. N. Spergel et al., Astrophys. J. 
Suppl. 170, (2007) 377. 

5. D. Clowe et al, astro-ph/0608407; D. Clowe, S. W. 
Randall and M. Markcvitch, astro-ph/0611496. 

6. G. Jungman, M. Kamionkowski, and K. Griest, Phys. 
Rep. 267, (1996) 195; G. Bertone, D. Hooper, and J. 
Silk, Phys. Rep. 405, (2005) 279. 

7. M. W. Goodman and E. Witten, Phys. Rev. D 31, 
(1985) 3059; I. Wasscrmann, Phys. Rev. D 33, (1986) 
2071; K. Griest, Phys. Rev. D 38, (1988) 2357; P. F. 
Smith and J. D. Lewin, Phys. Rep. 187, (1990) 203. 

8. A. K. Drukier, K. Preese, and D. N. Spergel, Phys. 
Rev. D 33, (1986) 3495; K. Freese, J. Frieman, and A. 
Gould, Phys. Rev. D 37, (1988) 3388. 

9. A. M. Green, JCAP 0708, (2007) 022. 

10. M. Drees and C. L. Shan, JCAP 0706, (2007) Oil. 

11. J. Engel, Phys. Lett. B 264, (1991) 114. 



